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Abstract 

We define an infinite set of families of graphs, which we call p-wheels and denote {Wh)^\ 
that generalize the wheel (p = 1) and biwheel (p = 2) graphs. The chromatic polynomial for 
(Wh)ffl is calculated, and remarkably simple properties of the chromatic zeros are found: 
(i) the real zeros occur at q = 0, 1, ...p + 1 for n — p even and q = 0, 1, ...p + 2 for n — p odd; 
and (ii) the complex zeros all lie, equally spaced, on the unit circle \q — (p + 1)| = 1 in the 
complex q plane. In the n — > oo limit, the zeros on this circle merge to form a boundary 
curve separating two regions where the limiting function W({(Wh)^},q) is analytic, viz., 
the exterior and interior of the above circle. Connections with statistical mechanics are 
noted. 
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1 Introduction 



The chromatic polynomial P(G, q) of a graph G plays a fundamental role in graph theory; 
it gives the number of ways of coloring this graph G with q colors such that no two adjacent 
vertices have the same color fl|]-0. Since P(G,q) < q n for a graph G with n vertices, it is 
natural to consider the asymptotic limiting function 

W({G},q)=\un P{G,q) l ' n (1) 

where {G} denotes the limit as n — > oo of the family of n-vertex graphs of type G. There 
is a natural generalization, which we assume here, of the variable q from integer to complex 
values. Recently, we calculated the W({G},q) functions for several families of graphs {G} 
H and studied their analytic properties. One of the interesting results that we found was 
that for three families of graphs, namely, circuit, wheel, and biwheel graphs the zeros of the 
respective chromatic polynomials, i.e., the chromatic zeros of these graphs, lie at certain real 
integer values and on certain unit circles in the complex q plane. A question that occurred 
to us was the following: Can one generalize this finding and if so, can one thereby achieve 
a deeper understanding of the result. We have succeeded in finding a comprehensive gener- 
alization which does yield further insight into this result, and we present this here. Given 
the elementary relation that the chromatic polynomial P(G, q) = Z(G, q,T — 0)paf, where 
Z(G,q,T = 0)paf is the partition function of the zero-temperature g-state Potts antiferro- 
magnet (AF) on the graph G, and the consequent relation, for the thermodynamic limit, 
between W ({G}, q) and the exponent of the ground state entropy of the Potts AF, our results 
also have implications for statistical mechanical properties of the Potts antiferromagnet on 
a certain family of graphs. 

To construct the p-wheel with n vertices, we start with a circuit with n—p vertices, C n _ p , 
and adjoin p points to each of the points on the circuit and to each other. Here, obviously, 
n — p > 1 and, for a nontrivial circuit, n — p > 2. This process could be symbolized as the 
product K p x C n _ p , where K p is the complete graph on p vertices. We denote the p-wheel as 
(Wh)^. This generalizes the definitions of the n-vertex wheel and biwheel graphs (Wh) n 
and U n ; these are seen now as special cases of the p- wheel for the respective values p — 1 
and p = 2: (Wh)^> = (Wh) n and (Wh)^ = U n . The n-vertex circuit graph C n is subsumed 
as the p = special case: (Wh)^ = C n . Thus, each family K p x C n - P = (Wh)^ forms an 
infinite set as one increases the number of vertices n. Clearly, 

(Wh) { ;l = K p+1 (2) 
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2 Theorems on Chromatic Zeros of (Wh) 



(p) 



n 



In general, the chromatic polynomial P(G, q) for a graph G has zeros at the discrete values 
q = 0, q = 1, and, if G contains one or more triangles, also the value q = 2. Furthermore, 
since P(G, q) is real for real q, it follows that the zeros of this polynomial (i.e. the chromatic 
zeros) are invariant under the complex conjugation q — > q*. We consider an n- vertex graph 
G in a particular family {G} of graphs for which, as n — > oo, aside from the above-mentioned 
discrete zeros at q = 0, 1 and, for G D A, q = 2, the remaining zeros merge to form the union 
of boundaries B separating various regions in the complex q plane. 0] A general question 
that one may ask about this latter set of zeros is whether all, or a subset of them, lie exactly 
on the boundary curves B even for finite n. (Some of the set of discrete zeros may also 
lie on B.) Clearly, the fact that a subset of zeros merge to form the boundary B in the 
n — > oo limit does not imply that, for finite graphs G, some subset of zeros will lie precisely 
on B. We have investigated this question and have found that there do exist some families 
of graphs {G} for which the zeros of P(G, q) (aside from a certain discrete set) lie exactly 
on the respective boundary curves B. Furthermore, we find that these families of graphs 
also provide an answer to the question posed at the beginning of this paper, concerning a 
generalization of the property that B comprises a (unit) circle. 

Our main results are contained in the following two theorems (as noted above, for a 
nontrivial circuit C n _ p in {Wh)^\ n — p > 2). 

Theorem 1 

The chromatic polynomial for the p-wheel is 

P{{Wh)i\ q) = [ l[(q -£)][(q- P - l) n ~ p + {-l) n - p {q - P - 1)] (3) 

1=0 



Proof 

We use the theorem that if a graph H is obtained from a graph G by adding a vertex which 
is adjacent to all of the vertices of G (i.e. connected to them by bonds of the graph), then 
(e.g., Ref. i) 

P(H,q)=qP(G,q-l) (4) 
We take G = C n _ p , the circuit with n — p vertices, and apply (f|) iteratively p times to get 

P((Wh)%\q)=[li(q-£)]P(C n - p ,q-p) (5) 

£=0 
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Substituting the basic result P(C n , q) = (q - l) n + (-l) n (q - 1) then yields (|). □. 
Next, we have 

Theorem 2 

The real zeros of the chromatic polynomial of the p-wheel, P((Wh)^\ q) occur at q — 
0, 1, ...,p + 1 for n — p even and q = 0,1, ...,p + 2 for n — p odd, while the complex zeros all 
lie on the unit circle \q — (p + 1)| = 1 in the complex q plane. Furthermore, the zeros are 
equally spaced around this circle, and in the n — > oo limit, the density g({(Wh)^}, 9) of 
zeros on the circle q = p + 1 + e t0 , —tt < 9 < tt, is a constant, independent of 9. 



Proof 



Defining x = q — (p + 1), we can write the second expression in eq. (fty as 



xix 71 ^- 1 + (-l) n - p ) 



(6) 



It follows that the complex zeros of P((Wh) , q) are the complex solutions to the equation 
x n-p-i _|_ ^_-Qn-p _ q ? ^ich clearly lie on the unit circle \x\ = 1, i.e., \q — (p + 1)| = 1. 
Specifically, for n—p even, x = —1, i.e., q = p is a root, and the n — p — 2 complex roots are 
given by 

g = p+l + e «-(p+i) , s = 1, .., n — p — 2 (7) 

while for n — p odd, x = ±1, i.e., q = p, p + 2 are real roots, while the n — p — 3 complex 
roots are 



g = p + 1 + e™-(p +1 ) 



s = 1, .., n — p — 2 ;s ^ 



n — p — 1 



where the (integer) value (n ~ p ~ 1> is omitted because this root is real. These complex roots 
obviously come in complex conjugate pairs; a factorization which makes this explicit is, for 
n — p even, 



P{{Wh)V,q) 



n—p even 



p+1 

11(9-4 

£=0 



— 2~~ ^ 



II |g-(p+l + e"-(p+ 1 ))||g-(p + l+e "-ip+d )| (9) 



-(2j+l)7Ti. 



and for n — p odd, 



P((W^)(f\g) n _ p odd 



p+2 

nc?- 

£=0 



Tl — p — 1 



tl {g-(p+l+e"-^+ 1 ))}{g-(p+l+e"-S+ 1 ))} (10) 



2jiri 
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The complex factors in eqs. (§) and (|10J) follow immediately from eqs. (0) and @. 
This proves the second part of our theorem, that the complex roots lie on the unit cir- 
cle \x\ = 1 with equal spacing. In particular, in the n — > oo limit, it follows that the density 
g({(Wh)^}, 9) of zeros on the circle q = p + 1 + e ie , —n < 9 < n is a constant, independent 
of 9. If one normalizes this density according to 

g({(Wh)M},9)d6 = l (11) 

J —IT 

then 

g({(Wh)^},9) = ±- (12) 

Concerning the product of real integer roots, we observe that for n — p even, eq. (^|) has 
the real factors x and x + 1, i.e., q — (p + 1) and q — p which, combined with the factor 
n£=o(<? — ■O i n ©> yield the factor n?=o(? — ■O i n (§)• F° r n — p odd, eq. (^j) has the real 
factors x, x + 1, and x — 1, i.e., g — (p + 1), g — p, and q — (p + 2) which, combined again with 
the factor n£j(g - £) in (|), yield the factor n?=o(? - in ©■ This proves the first part 
of the theorem, which specifies the real roots of P((Wh)^\ q). This completes the proof of 
both parts of the theorem. □ 

3 Analytic Structure of W({(Wh)®}, q) 

We next consider the limiting function W({(Wh)^},q). As we have discussed in Ref. 0, 
eq. ([!]) is not, in general, adequate to define this function because of a noncommutativity 
that occurs at certain special values of q, which we shall denote q s , where 

lim lim P(G,q) 1/n ^ lim lim P(G,q) 1/n (13) 

n— >oo q— *q s q—"ls n—^oo 

The origin of this noncommutativity of limits is an abrupt change in the behavior of P(G, q) 
in the vicinity of such a point q s ; for q ^ q s , P(G,q) grows exponentially as the number 
of vertices n in G goes to infinity: P(G, q) ~ a n for some nonzero a, whereas precisely at 
q = q s , it has a completely different type of behavior; typically, P(G, q s ) = co(q s ) where Co(q) 
may be a constant, independent of n or may depend on n, but in a way that does not involve 
exponential growth, like (—1)™. The set of special points {q s } includes q — 0, q — 1, and, if 
G D A, also q = 2; at these points, P(G, q s ) = 0. 

In the present case, we find that for a fixed p, in the n — ► oo limit, the zeros on the 



unit circle \q — (p + 1)| = 1 merge to form a (complete) circle with a density given by (12). 
This circle thus constitutes the boundary curve B and divides the complex q plane into two 
regions: 

Ri : q such that \q-(p + 1)| > 1 (14) 
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and 

R 2 : q such that \q — (p + 1)| < 1 (15) 

The set of special points {q s } associated with the noncommutativity in ( |T3D are the real 
integral zeros of P((Wh)^\ q), i.e., {q s } = 0, 1, ...,p+l for n—p even and {q s } = 0, 1, ...,p+2 
for n — p odd. We find that 

W({(Wh) ip) }, q) = q-(p + l) for q G R x and g £ {g s } (16) 

Because of the noncommutativity (0), the formal definition ([I]) is, in general, insufficient to 
define q) at the set of special points {q s }', at these points, one must also specify the 

order of the limits in flT3|). One can maintain the analyticity of ^({G}, q) at these special 
points q s of P{G, q) by choosing the order of limits in the right-hand side of eq. (|TB|): 

W({G},q s ) Dqn = lim limP(G,g) 1 /" (17) 

(J — '(/s ll — 'OO 

As indicated, we shall denote this definition as D qn Although this definition maintains the 
analyticity of W({G},q) at the special points q s , it produces a function W({G},q) whose 
values at the points q s differ significantly from the values which one would get for P(G, g s ) 1//n 
with finite-n graphs G. The definition based on the opposite order of limits, 

W({G},q s ) Dnq = lim lim P(G,g) 1 /" (18) 

gives the expected results like W(G, q s ) = for q s = 0, 1, and, for G D A, q = 2, but yields 
a function W({G},q) with discontinuities at the set of points {q s }- In the present case, we 
have 

W({(Wh)^},q) Dnq =0 for q E R 1 and q G {q s } (19) 

while 

W({(Wh)^},q) Dqn =q-(p + l) for q G R 1 and q G {q s } (20) 

For g G i?2 there is no canonical choice of phase to use in the n'th root ([[]) |J; however, for 
the complex amplitude we have 

\W({(Wh) ip) },q)\ = 1 for g G P 2 and g £ {g s } (21) 

The region i?2 contains the point g = p + 1 which is a member of the set of discrete 
zeros {q s }; for g = p + 1, the two orders of limits give W({(Wh)^},q)D nq = and 
\W({(Why^},q) Dqn \ = 1. The limiting function W({(Wh)^},q) is nonanalytic on the 
boundary circle \q — (p + 1)| = 1. 
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In Ref. || (see also Ref. ||) we discussed the quantity q c ({G}), the maximal real point 
of nonanalyticity of q), and inferred the values of q c for certain lattices {G} = A. A 

corollary of our theorems above is that 

qc({(Wh)M})=p + 2 (22) 

As was the case with the lattices previously considered, we infer that for integer q above 
this critical value, i.e, here, q > p + 2, the Potts antiferromagnet on the infmite-n limit of 
the p-wheel graph has a nonzero ground state entropy Sq = In W({(Wh)^}). For integer 
q < p — 1, the model is frustrated because of the impossibility of assigning different colors 
(i.e., spin states in the Potts AF) to all of the adjacent vertices of the K p subgraph. It 
would be interesting to investigate further the properties of the Potts antiferromagnet on 
the infinite-n limits of p-wheel graphs. 

Because of the connection to statistical mechanics, it is worthwhile to contrast our theo- 
rem with the famous Yang-Lee circle theorem || . Yang and Lee considered the Ising model 
partition function Z = J2{a r } e~^ n at temperature T in an external magnetic field H, for 
which the Hamiltonian is Ji = — J2 r ,r' <J r Jr,r'&r> — HJ2 r °V, where oy = ±1 and (3 = (/c^T) -1 . 
They proved that for an arbitrary (finite as well as infinite) lattice with spin-spin couplings 
which are ferromagnetic, J r y > 0, but not necessarily nearest-neighbor (and physical tem- 
perature < T < oo), the zeros of Z in the variable /i = e~ 2/3H lie on the unit circle \p\ = 1 
in the fi plane. In the thermodynamic limit, these zeros merge, becoming dense and forming 
an arc fi = e l6 for ^o(^) — \@\ < 7r, where ±#o(T) denote the angles at the endpoints of this 
arc. The free energy is nonanalytic across the arc. For temperature T > T c , the arc leaves 
a gap which includes the positive real axis in the fi plane, but as T decreases through the 
critical temperature T c , this arc closes, pinching the real axis at fi = 1. For < T < T c 
the zeros continue to form a full closed unit circle. Our theorem shares with the Yang-Lee 
theorem the property of zeros lying on a unit circle and of being true for finite graphs as 
well as in the limit n — > oo. However, there are several important differences: (i) first, in 
the Yang-Lee circle theorem, all of the zeros of Z lie on the fixed unit circle = 1, whereas 
in contrast, in our theorem, only a subset do and, moreover, the center of the circle is at 
a different position, q = p + 1 for each respective family of p-wheel graphs, (Wh)^; (ii) 
the Yang-Lee theorem describes a multivariable problem, dealing with the zeros of Z in the 
complex fi plane, which depend parametrically on the temperature T, and it requires in its 
premise that all J r y > 0, whereas in our case, once p is specified, i.e., one specializes to 
a given p-wheel graph (Wh)^\ the zeros in the complex variable, q are fixed; (iii) thus, 
while in the n — > oo limit, the Yang-Lee zeros merge to form an arc which does not form 
a closed circle for T > T c , the subset of chromatic zeros that we show lie on the unit circle 
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\q — {jp + 1)| = 1 always close to form a complete circle in the n — > oo limit; and (iv) on the 
Yang-Lee arc, the density of zeros is not, in general, a constant, independent of 9, whereas 
the density g{{{Wh)^}, 9) in our case is, in fact, independent of 9. Despite these differences, 
the properties that our theorem and the Yang-Lee circle theorem share in common, namely 
that they both apply to a finite graph and imply zeros lying on a unit circle, show a very 
interesting connection between these results. 

4 Conclusions 

In conclusion, we have constructed a set of families of graphs which we call p-wheels, {Wh)^\ 
and calculated their chromatic polynomials. We have proved that these provide a (countably 
infinite) generalization of the behavior that we had found previously for the circuit, wheel, 
and biwheel graphs (which are seen as the special cases p = 0, 1 and 2, respectively), namely 
that a subset of the chromatic zeros lie exactly on the unit circle \q — (p + 1)| = 1 in the 
complex q plane, and in the limit as the number of vertices n — > oo, these zeros merge to 
form the full circle \q — (p + 1)| = 1 with density independent of position on the circle. Since 
the chromatic zeros lie exactly on this circle even for finite n, the families also constitute a 
constructive answer to another question that we posed, namely to find families of graphs such 
that chromatic zeros lie precisely on the asymptotic boundary curve B. We have compared 
our results with the Yang-Lee circle theorem and remarked on the connections with zero- 
temperature Potts antiferromagnets. Our present theorems provide a deeper and unified 
understanding of our previous findings in Ref . . 

This research was supported in part by the NSF grant PHY-93-09888. 
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